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Abstract 

The ground state energy of ideal a-matter at T = is analyzed within the 
framework of variational theory of Bose quantum liquids. Calculations are 
done for three local a — a potentials with positive volume integrals and two- 
body correlation functions obtained from the Pandharipande-Bethe equation. 
The energy per particle of a matter is evaluated in the cluster expansion 
formalism up to four-body diagrams, and using the HNC/0 and HNC/4 ap- 
proximation for a Bose liquid. At low densities the two methods predict 
similar EOS whereas at higher densities they are sensitively different, the 
HNC approximation providing saturation at lower density, bellow the satu- 
ration value of nuclear matter. Inclusion of higher order terms in the cluster 
expansion of the condensate fraction is leading to a stronger depletion of the 
alpha condensate with the density compared to the two-body approximation 
prediction. 
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1. Introduction 

Renewed interest in the properties of a matter is manifest in the literature 
especially in connection with a-particle Bose-Einstein condensation (BEC) 
in a- like nuclei (see [H and references therein). Calculations reported in 
this reference are pointing to the existence of a Bose-Einstein condensate of 
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a-particles at low densities. It was also noted by these authors that with 
increasing density the condensate fraction is reduced such that at density 
corresponding to the saturation of nuclear matter 0.04 a particles per 
fm^), the condensate fraction is reduced to roughly one half. The estimation 
of the condensate fraction was done in the lowest approximation, i.e. the 
radial distribution function (RDF) is approximated by the square of the 
two-body correlation function (CFN), and therefore it is less justified for 
higher densities. There is howewer an old estimation by Clark and Johnson 
[2I for three values around the saturation density of nuclear matter using the 
hypernetted chain approximation in the lowest order (HNC/0), i.e. taking 
into acount only nodal diagrams in the infinite density expansion of the RDF. 
It provides a severe reduction of the condensate fraction (;^15%) compared 
to the lowest-order cluster expansion at the same density. On the other 
hand calculations of the cold a matter equation of state (EOS) reported 
by the same authors within the HNC/0 approximation and using the soft 
core a — a potential of Ali and Bodmer [3| are predicting the saturation 
point at a high density (p^ ~ 0.085 a particles per fm^). These benchmark 
calculations of the a matter EOS were very recently compared to results 
obtained in the frame of the scalar 0^ effective field theory with negative 
quartic and positive sextic interactions, to simulate the attractive character 
at long distances and repulsive at short distances, and found to be in a very 
good agreement Though no estimations of the condensate fraction are 
provided for this high density saturation point, from the estimation made 
at lower densities, as quoted above, we expect a stronger depletion of the 
BEC. One is then confronted with the problem that at low densities the a 
matter condensate is far from equilibrum, whereas at the saturation point 
the condensate fraction is small. 

It was advocated that beyond a critical density (p* ~0.03 nucleons per 
fm^), due to the strong overlap of the wave-functions and the unavoidable 
action of the Pauh principle, a total extinction of the a structure should occur 
The phenomenological a— a potentials used in the past are systematically 
predicting saturation of a matter at densities considerably larger than this 
critical densities. It would then be important to establish if the saturation 
of the a matter takes place below this critical density if one employs other 
types of potentials that incorporate more microscopic input. 

The aim of this letter is to analyse the a matter EOS over a wide range of 
densities and try to find the optimal CFN which reflects the interplay between 
the strong short-range and the long-range correlations that ultimately would 
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lead to saturation. Gaining insight in the saturation properties of a matter 
could also shed hght on the condensate fraction reduction issue. 

In what follows the g.s. energy of an infinite system of neutral a-particles 
interacting via two-body forces is calculated within the variational theory of 
Bose liquids. As input to the energy calculation we use a prescriptions for the 
CFN obtained by extremizing the energy functional in the two-body cluster 
approximation. The g.s. energy is then calculated via the cluster expansion 
adding the three-body and four-body correlations and with the HNC method 
that is more reliable in the high density sector. 

2. Two-Body Potential 

The a — a potential is necessarily nonlocal due to the short range repul- 
sion between a particles coming from exclusion effects. Phenomenological 
potentials obtained by inversion suffer of some difficulties which have not 
been yet fully resolved, mainly due to the limited range of energies where 
phase shifts have been measured. We remind the reader that potentials ob- 
tained by inversion are unique if and only if the corresponding phase shifts 
are known for fixed angular momentum for all energies, up to infinity, as 
required by Marchenko theorem j^. If there exist bound states, asymptotic 
normalization coefficients for these states are necessary. Further, a — a po- 
tentials are required to reproduce at least qualitatively the known resonances 
in *Be. These difficulties lead to a lot of ambiguities in the proposed a — a 
potentials. For the typical Ali-Bodmer potential the nonlocality translates 
into a strong dependence of angular momentum in order to reproduce the 
repulsive effect of the redundant states in the L=0 and L=2 partial waves. 
We shall use the S-state Ali-Bodmer potential as representative for the class 
of phenomenological potentials with a soft inner repulsive core and a weak 
long range atractive component. There is however a more stringent require- 
ment to be imposed on a — a potential for a matter calculations. In ref. 0] 
the Euler-Lagrange equation for a Bose system within the HNC approxima- 
tion was analized in detail and shown that a spherically symmetric solution 
which do not lead to a colapse of the system requires a potential with positive 
volume integral. 

The Ali-Bodmer potential the a matter actually fails to saturate. In 
fact a very shallow minimum in the EOS at a high density is predicted. 
Other schematic potentials (hard core) are trivially saturating at densities 
and energies close to the nuclear matter saturation point {pa ~ 0.04 a/fm^. 
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and E/Na ~ -11-16 MeV). However as we noted above, at such high densities 
the a-condensate is almost completely depleted according to the variational 
approach. Somehow this dissapointing result is conflicting with what would 
one expect based on the manifestation of a clustering in light real nuclei. 
The clusterization of a particles on the surface of nuclei at densities around 
half the central nuclear density, as revealed by a-decay, a-transfer reactions 
or the putative dilute three-alphas condensate in the Hoyle state of ^^C are 
pointing to a higher stability of a matter at lower densities. 

Since the potentials providing saturation at lower densities are highly 
schematic (infinite repulsive short-range interactions) we resort to a calcula- 
tion of the bare a — a interaction based on the double-folding method for two 
ions at energies around the barrier. As input we consider realistic densities 
of the a-particle and modern effective nucleon-nucleon interactions. 

In the double-folding framework [sl the interaction between two alpha 
clusters is calculated as a convolution of a local two-body potential Vnn and 
the single particle densities of the two clusters 

Vaa{r) = j dridr2pa{ri)pa{r2)vnn{p, r -ri + r-2) (1) 

The effective n — n interaction Vnn is taken to depend on the density p 
of the nuclear matter where the two nucleons are embedded. It should also 
consist of a density independent finite-range part with preferably two ranges 
such that a potential similar to the Ali-Bodmer is obtained. We therefore 
choose the Gogny effective interaction. Out of the three main parametriza- 
tions of this interaction, only the Dl [ij and the most recent one DIN l3 



are satisfying the requirement of positive volume integral of the a — a poten- 



tial. The DIS parametrization [ll[ leads to a potental with negative volume 
integral and is omitted. . 

In what follows we include only the direct part in the double-folding po- 
tential ([1]). The knock-on nonlocal exchange component leads to a strongly 
attractive local component mainly due to the fact that the Perey-Saxon lo- 
calization procedure is not reliable at the energies of interest considered in 
this paper. 



We take a Gaussian nuclear matter distribution inside the a-particle 
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P„(r)=4(— e-^V^^ (2) 
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with an oscillator parameter b that corresponds to a root mean square (rms) 
1.58±0.004 fm resulting from a Glauber analysis of experimental interaction 
cross sections 
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The direct effective n — n force in the Gogny parametrization reads : 



3 

+ 



i=l 

\{rr~r2) (3) 



Inserting the gaussian density distribution ([2]) in the double folding integral 
([1]) and using the Campi- Sprung prescription [13| for the overlap density 

P(l, 2) = (^p„(ri - is)p,(r2 + ^s)^ ' , (4) 

where s = ri + r — r2 is the n — n separation in the heavy-ion coordinate 
system jsf, we obtain the a — a potential, 

2 / 2 \3/2 

v.^{r) = aY^{AW, + 2B.,-2H.,-M;)^-^^^ e~^V,l+2i>- 

3 4'^+^ _7+2 2 

+ -2'j^i^W^K7W^'~^^' 

which is assumed to be the same in all partial waves. In the left panel of Fig{T] 
we represent the three potentials on a magnified scale around the minimum. 
The two selected Gogny interactions display shalower pockets that are shifted 
to larger radii compared to the Ali-Bodmer potential. Before ending this 
section we remark that a fictions ^Be (Coulomb interaction is switched off) is 
slightly bounded by Ali-Bodmer {J^ = 0+, £'=-1.63 MeV), loosely bounded 
by Gogny-DIN (J^ = 0+, E=-Q.16 MeV) and not bounded by Gogny Dl 
interaction, in agreement with the patterns depicted in Fig. [T] Levinson's 
theorem applied to the a — a system [14] which predicts So{E = 0) = Sn, 
where 6o is the scattering phase for the L = channel, is reasonably well 
satisfied by all potentials. 

3. Variational approach 

Within the variational approach to the description of a uniform system 
of spinless Bose particles at zero temperature interacting via a two-body 
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potential v, 



the expectation value of the energy 



E\ 





H 









(6) 



can be calculated using the Bijl-Dingle-Jastrow (BDJ) trial wave function 

Above, /(r) is a positive semi definite function that is required to reflect 
exclusively the two-body correlations. Thus at short distances, where the 
potential is expected to be strongly repulsive, / is small , whereas at large 
distances (r — > oo), when the two-particle " decorrelate" , /(r) — ^ 1. 

The Jackson- Feenberg energy £^jf, measured relative to the rest energy of 
Na, a-particles, is [ly, [l7| 



£ 



JF 



j drg{r) 



2m, 



-A In /(r) 



(8) 



where g{r) is the two-body radial distribution function (RDF). This form of 
the energy has the advantage for bosons that terms including the three-body 
radial distribution function does not occur. 
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The cluster expansion formalism heavily relies on the fact that h = p — 1 
is of short-range and consequently the integral uj = J h{r)dr is small com- 
pared to the volume Vt of the system (a; = —\ = « [isj. Then g can 



be evaluated by an expansion in powers of p [19||. Including up to four-body 
diagrams the expansion reads 

g{ru) = f (ri2)|l + py"t/r3/i(ri3)/i(r23) 

+ ^p^ j drs j dri[2h{rri)h{r2i)h{ru){l + 2h{ri^)) 

+ /i(ri3)/i(r23)/i(ri4)/i(r24)(l + M^34)] + 0(p')} (9) 

The last term in the above formula include ring, diagonal, open and elemen- 
tary four-body diagrams. 

Retaining only the first term in the cluster expansion of g, the lowest-order 
Jackson-Feenberg functional is obtained upon substitution in the expression 
dHD of £ 

^2[/, V/, A/] = \pjdr ^vUr)f{r) + ^ [{yf{r)f - /(r)A/(r)] 

(10) 

A practical method to obtain / is provided by the Pandharipande-Bethe 
prescription 20|] . It consists in varying the two-body energy functional with 
respect to the one-parameter family of functions / under the constraint of 
normalization p j[p — l\dr = —1. Thus, the variational problem is reduced 
to the eigenvalue problem 

(d'f 2df\ 

[-r^ + -^r] + {Va^o.- X) f = Q 11 

iTLa \ dr"^ r dr J 

The CFN is subjected to natural boundary conditions that reflect the 
Schrodinger condition at the origin of a soft-core and the necessity to heal 
to unity at a given distance d 

u{r) ^ rf{r) ^ 0, /(r > rf) = 1, f{d) = (12) 

For each density we varried the parameter d untill the normalization condi- 
tion was fuUfiled. In this way we obtained a solution that is density dependent 
(since we take d ~ p'^/'^), and the overshoot (the peak exceeding the unity) 
increase with density as can be seen in Fig|2]for all three potentials. 
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Figure 2: Correlation function f{r) resulting from solving eg. pTjl for the two cases : a) 
dependent on density at three different densities and b) unique for all densities. Each 
figure contains these correlations functions for a given a ~ a potential. 

Another recipe that we apphed was to relax the normahzation constraint 
and match instead the potential at the right boundary with the eigenvalue 
A = Vaa{d)- This constraint is fullfilled for a unique d regardless the value 
assigned to p. The interpretation assigned by Pandharipande to the Lagrange 
multiplier A is that of a contribution to the average field comming from the 
excluded particles. The corresponding CFN is then density independent 
and it does not display any overshoot. In both cases we select the lowest 
eigenvalue, corresponding to zero nodes of the CFN. 

Naturally the validity of the cluster expansion is limited to low densities 
if the truncation is performed after including the four-body correlations. A 
method that allows the determination of g and is more reliable for higher 
densities is provided by the HNC approximation, where the diagrams be- 
longing to the nodal set A/" are summed up to infinity. Thus, the RDF in 
HNC/0 approximation is defined as j2l[ 

g{r) = /2(r)e^W (13) 

The function N{r) results from the iterative solution of a set of non-linear 
equations 

X(r) = ^(r)-A^(r)-l (14) 
N{r) = j dridr2X{ri){X{r2) + N{r2))5{ri + r -r2) 

(15) 
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Within the HNC/0 approximation elementary diagrams are neglected. 
The first correction to this approximation is given by HNC/4 , where the 
label 4 indicate that only the four-body elementary diagram is taken into 
account 



22|. As a consequence expression f[T3|) is modified, such that this 



diagram is added to the infinite number of nodal diagrams : 

g(r) = /(r)2e^M+^*M (16) 



where 23 



e4(ri2) = ^p' j rfr3rfr4G(ri3)G(r23)G(ri4)G(r24)G(r34) (17) 

and G{r) is a short-hand notation for g{r) — 1. Thus, within HNC/4 approx- 
imation we obtain a fully selfconsistent solution of equations (fT^ - (fT7|) . 



4. Ground-state energy of ol matter 

The cluster expansion of the g.s energy results from inserting the expan- 
sion in powers of densities of the RDF feq.dp])). The two-body and three-body 
terms are already given in the literature [19|. The four-body terms can be 
easily worked out employing the folding technique. Explicit expressions of 
the four independent diagrams contributing to this term were derived in . 
On the other hand the g.s. energy in the HNC/0 approximation is obtained 
by simply substituting into ([H]) the RDF obtained from the iterative solution 
of eqs.dnD-dTS]) 

Fig. 3 displays the energies for the three potentials Ali-Bodmer, Gogny Dl 
and DIN calculated within the cluster expansion method (solid curves) and 
with the HNC/0 method (circled crosses) for both CFN resulting from the 
Pandharipande-Bethe equation. Upon comparison with the a — a potential 
we remark a similar pattern. It became obvious earlier that when we move 
from the most attractive potential (Ali-Bodmer) to the most repulsive one 
(Dl), the "molecular" pocket is shifted towards larger distances and shal- 
lower minima. The saturation point of the a matter EOS is considerably 
displaced to lower densities and instead of the overbinding manifest for the 
Ali-Bodmer potential we find a low-binding for the Gogny-Dl case. Another 
important feature of the obtained EOS is that the HNC/0 calculations pro- 
duce less bound a matter with an obvious preference to saturate at lower 
densities. Although the cluster expansion formalism predicts softer EOS in 
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Figure 3: Total energy calculated with the cluster expansion (CE) taking into account up 
to four-body correlations contributions (solid lines) and with the HNC /O method (circles) 
for the " /9-independent" CFN considered in this paper. 



the case of the p-independent solution, the HNC/0 resuhs are almost iden- 
tical for the two different choices of the CFN, a fact that suggests that most 
likely both are close to the optimal solution. As noted in 2J| for the Bose 
homework problem, the fact that different looking correlation functions may 
give extremely close energies, is a general feature of the variational method. 

The cluster expansion and HNC methods are similar over a larger range 
for the Ali-Bodmer potential, i.e. up to p ~ 0.0225 a per fm^ whereas for 
the Gogny-Dl potential the interval is limited to p < 0.0125 a per fm^. 

Previously HNC/4 calculations have been reported for the ^He quantum 



liquid [23|, |25| but to our knowledge there is no similar calculation for the 
g.s. of a matter. In FigJUwe compare the EOS predicted in the two HNC 
approximations using the Gogny potentials. Apparently the contribution of 
the 64 diagram results only in a small lowering of the saturation energy. 

From the EOS derived within the HNC/0 formalism we extracted the cor- 
responding incompressibilities as a function of density (see FigHj). Although 
the EOS for the Gogny-Dl interaction has the largest curvature, since the 
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Figure 4: Comparison between the EOS predicted with HNC/0 and HNC/4 for Dl and 
DIN potentials for both CFN selections. 



saturation takes place at a rather small density (psat ~ 0.015 a per fm^), 
the incompressibility is also small , 64 MeV. Instead the Gogny-DIN 

interaction which predicts a saturation point at approximately psat 

^ 0.0325 

a per fm^ and binding energy -16.8 MeV, provided we add the -7.07 MeV 
binding energy per nucleon of the a particle, has a corresponding value of 
K close to the generally accepted range of standard nuclear matter values, 
i.e. K fti 265.4 MeV. We note that within the dynamical lattice model of 
a matter with a Brink-Boeker I effective interaction by Tohsaki, a volume 
energy of -16.76 MeV and incompressibility K ^ 181 MeV are predicted 26 . 

From the inspection of FigJUwe conclude that the interval where a matter 
develops thermodynamical instability {K < 0) is ranging over dilute densities 
for the Gogny Dl and DIN potentials, whereas for the Ali-Bodmer interaction 
also large values of the densities are affected by instability. 

In our opinion the domain of thermodynamical instability provides a fur- 
ther ground to discard the Ali-Bodmer potential since it predicts a highly 
unstable a matter at densities where a clusterization is believed to be im- 
portant. The Gogny Dl potential has saturation properties very close to the 
symmetric nuclear matter as can be noticed on FigHl On the other hand 
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although the Gogny-Dl interaction predicts a much softer EOS and a low 
binding energy, it displays instead stability at densities not far from the Mott 
density We therefore expect a realistic a matter EOS to develope a satu- 
ration point approximately in the domain bounded by the saturation points 
of the Dl and DIN EOS. 




Figure 5: Incompressibility of a matter within HNC/0 for the three interactions Ah- 
Bodmer, Dl and DIN. The fuU circles specifies the saturation points and the dashed 
region the accepted range of symmetric nuclear matter saturation values. 



5. Reduction of the Condensate Fraction 

Another quantity of interest for the a matter ground state is the fraction 
of particles nc filling the zero momentum state (BEC) . Using a diagramatic 
expansion, Ristig and Clark showed that nc can be represented by irreducible 
n-body cluster diagrams AQ^"^ with two external points and n — 2 internal 
field points [23] : 

no = e^«'^'+^«'^'+^Q<^'+- (18) 

For the first three terms closed forms are available. Introducing the notation 
({r) = /(r) — 1 and making use again of the compact representation of 
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integrals via the convolution operation, the corresponding expressions are 



Ag(i) = pj drC{r) (19) 



AQ(2) =p^ j dr |c(r)(C * G){r) - hi{r){h * G)(r) 



(20) 



Ag(') = ^pM rfr2dr3rfr4(2C(ri2)C(n3)C(n4)-/^(ri2)/^(ri3)/^(ri4)) 



6' 

X [^3(^2, r-3, r,) - G(r23) - G(r43) - G(r24) - 1] (21) 

The calculation of the third-order contribution is greatly simplified by the 
fact that at the HNC/0 level, the tree-particle radial distribution function is 
given by the Kirkwood superposition approximation 

g^{r2, r3, r^) = g {r 23) g{r 43) g{r2i) (22) 

Note that the three-body RDF is large only when all three alphas are close 
to each other (within distances of order p~^/^). Details on the calculation of 
this last term will pe presented elsewhere. 

The exponential reduction of the condensate fraction is represented in 
Fig. [5] for the case when the second-order and the second plus third-order 
terms are added to the lowest order contribution for the three potentials 
used in this work. We perform an exemplification using only the density- 
dependent CFN. In all cases the third-order diagram is only slightly reducing 
Tic whereas the contribution of AQ^"^^ is important beyond a certain density 
that is approximately 0.01 a/fm^. As expected, in the case of the Ali-Bodmer 
potential the depletion of the condensate state with increasing density is 
smaller that in the case of Gogny potentials. Thus, whereas in the Ali- 
Bodmer case at saturation densities of nuclear matter 0.04 a/fm^) the 
depletion is around 35%, for the repulsive Gogny-Dl potential only a tenth 
of alphas are left in the condensate. On the other hand the HNC calculations 
employing Gogny potentials are predicting saturation at densities p ^0.015 
a per fm^ (Dl) and p ^0.0325 a per fm^ (DIN). As can be inferred from 
FiglSlthe condensate at this density is still well populated, i.e. 55% (Dl) and 
30% (DIN). 
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p (a particles per fm'"^) p (ci particles per fm"'') p {a particles per fni"^) 

Figure 6: Reduction of the condensate fraction due to the contribution of higher order 
terms in the cluster expansion of the function Q for the three potentials used in this work. 

6. Conclusions and Outlook 

The properties of cold a matter have been analyzed using the frame- 
work of the variational theory. Compared to previous studies, performed 
mainly by Clark and collaborators, we introduce two new a — a potentials 
that are based on the realistic a-particle densities and effective Gogny n — n 
interactions widely used in modern nuclear structure studies. Hypernetted 
chain calculations including the elementary 4-body diagram are providing 
slight corrections to the HNC/0 approximation. For the low-density regime 
the cluster expansion method is providing results similar to the HNC, which 
means that predicitions of the a matter g.s. properties in this case can be 
safely and rapidly obtained taking into account at most four-body correla- 
tions. 

To our knowledge for the first time in the literature a detailed investiga- 
tion of the condensate fraction for a matter was performed. The reduction 
of condensate fraction is enhanced once we take into account higher-order 
contributions in the cluster expansion of Q. For the Ali-Bodmer case the 
reduction is temperated by the more attractive character of the potential. 

We inferred that only the Gogny-Dl based potential is predicting a sat- 
uration of alpha matter close to the Mott density. Moreover at this density 
there is a significant condensate fraction. On the other hand the Gogny-DIN 
EOS displays saturation properties (density, energy and incompressibility) 
similar to the symmetric nuclear matter EOS. For this reason we are con- 
fident that a realistic saturation point of the a matter should be located 
approximately in the range between 0.015 - 0.0325 a per fm^. 

As mentioned in the introductory section the main factor that triggered 
this work was the renewed interest in the literature regarding the properties 
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of the a-condensate. We showed that compared to the dilute approximation 
adopted in ref. [l|, the inclusion of higher order terms in the cluster expansion, 
produce a steeper depletion of the condensate, which is even more accentu- 
ated for the repulsive Gogny potentials. 

The present work was also motivated by modern theories of low-density 
nuclear matter composed of neutrons, protons and alpha particles near the 
neutrinosphere during supernovae core-collapse (see |28|). In this context it 
would of interest to extend the present calculations to the case of non-zero 
temperature and inhomogenous distribution of a matter. 
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